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摘 要 : 针对 奇异 值 分 解 水 印 算法 导致 唐人 警 率 高 、 稳 健 性 不 强 的 问题 ， 提 出 一 种 基于 A A 4E E27 RE (NMF) fe 
强 奇 异 值 分 解 (BN-SVD) 相 结合 的 自 适 应 零 水 印 算 法 。 首先 将 原始 灰 度 图 像 进行 二 级 离散 小 波 变 换 (DWT), 对 变换 后 
的 二 级 低频 子 带 (LL2) 进 行 不 重 二 分 块 ， A AE E A 
采用 增强 奇异 值 分 解 , 依据 每 一 个 块 矩 阵 的 最 大 奇异 值 与 整体 最 大 奇异 值 均值 的 大 小 关系 构成 特征 向 量 ;利用 生成 的 
特征 向 量 与 经 过 Arnold 变 
自 适应 确定 增强 奇异 值 分 解 中 最 抗 攻击 缩放 比例 的 参数 有 。 实 验 结果 表明 , 在 虚 警 问题 上 NC 值 达 到 0.4 以 下 , JPEG 
pe 唆 声 、 滤 波 、 旋 转 、 剪 切 以 及 混合 攻击 下 ， 提 取水 印 图 像 与 原水 印 图 像 的 归 一 化 系数 NC 值 均 可 达到 99% vA 
上 ， 该 方案 高 效 地 解决 虚 警 问题 ， 具 有 和 较 强 的 稳健 性 ， 能 够 有 效 地 抵抗 各 种 攻击 。 
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Adaptive zero-watermarking algorithm based on block NMF and 
boost normed singular value decomposition 


Xiao Zhenjiu, Ning Qiuying, Zhang Han, Tang Xiaoliang, Chen Hong 
(College of Software Liaoning Technical University, Huludao Liaoning 125000, China) 


Abstract: For the problem of the singular value decomposition (SVD) watermarking algorithm caused high false alarm rate 
and low robustness, this paper proposed an adaptive zero watermarking algorithm , which was based on block NMF and 
boost normed singular value decomposition. First, the original gray image was transformed into two levels of DWT 
transform, then after the transformation of the LL2 subbands were divided into non-overlapping blocks, and each sub-block 
was broken down into NMF with rank r, then boost normed singular value decomposition was used for eigenmatrix derived 
from NMF decomposition. According to the relationship between the maximum singular value of each block matrix and the 
mean value of the global maximum singular value, it constructed feature-vector . Finally the generated feature vector was 
xor operation with the random encrypted watermark image of Arnold and the logistic map to generate a zero watermark. The 
parameter p in the singular value matrix must be determined adaptively by the optimization algorithm (BAS) , and find out 
the scaling ratio of the most resistant attack. The experimental results show that, in the false alarm problem, the NC value is 
below 0.4. Under the condition of JPEG compression, noise, filtering, rotation, shearing and mixed attack, the normalized 
coefficient NC of the extracted watermark image and the original watermark image can reach more than 99%. The method 
can solve the false positive errors efficiently and has strong robustness. Can resist all kinds of attacks effectively. 
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引言 域 水 印 算法 。 空 间 域 水 印 通常 是 直接 修改 图 像 的 灰 度 值 来 实 

m 山 效 果 。 变 换 域 水 印 则 是 通过 修改 载体 图 像 频 域 系数 来 达到 

近 20 年 来 ,数据 压缩 性 能 的 提高 和 互联 网 带宽 的 不 断 增 ”嵌入 目的 。 两 者 相 比 ， 变 换 域 水 印 更 鲁 棒 。 小 波 域 水印 大 多 
Jn. 使 得 数字 图 像 、 视 频 和 音频 更 容易 被 播放 , 复制 与 窜改 。 将 水 印 嵌 在 中 频 域 。 文 献 [3] 提 出 一 种 基于 小 波 变 换 和 NMF 
越 来 越 多 的 人 想 要 分 享 他 们 赋 有 创造 性 的 数字 作品 。 因 此 ， 的 数字 图 像 水 印 算法 ， 原 始 图 像 进行 二 级 离散 小 波 变 换 ， 选 
所 有 数字 作品 的 版 权 保护 和 所 有 权 申 报 的 都 变 得 非常 重要 。 择 图 像 中 频 区 域 进行 嵌 入 ， 对 LH2 做 NMF 分 解 ， 将 经 过 
为 了 保护 数字 作品 免 受 盗版 ， 数 字 水 印 提供 了 一 种 潜在 的 解 Arnold 变换 后 的 水 印 以 嵌入 强度 a 崔 在 到 分 解 后 的 特征 矩阵 
决 方案 。 因 其 在 数字 版 权 管理 和 保护 中 的 重要 应 用 ， 在 很 大 He 文献 [和 在 文献 [3] 的 算法 上 进行 了 改进 , 将 宿主 图 像 进行 
程度 上 推动 了 数字 水 印 研究 的 不 断 深入 。 三 级 小 波 变换 ， 选 择 低频 近似 分 量 进 行 嵌 入 ， 最 后 通过 量化 
数字 水 印 !3 是 指 将 信息 (水 印 ) 潜 移 默 化 地 嵌入 到 数字 文 ”系数 矩阵 得 到 表示 原始 图 像 的 特征 向 量 ， 然 后 通过 计算 水 印 
档 中 以 提供 内 容 保护 或 内 容 认证 的 过 程 。 按 照 嵌 入 式 算 法 分 ”和 特征 向 量 嵌 入 版 权 信息 。 引 入 离散 小 波 变换 后 的 算法 ， 水 
类 ， 主 要 可 分 为 两 种 : 一 种 是 空间 域 水 印 算法 ， 一 种 是 变换 印 的 不 可 见 性 得 到 提高 。 对 一 些 图 像 退 化 过 程 和 JPEG 压缩 
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的 鲁 棒 性 有 所 提升 。 因 单一 的 小 波 变 化 无 法 使 水 印 获得 强 鲁 性 ， 更 能 体现 智能 数据 处 理 的 本 质 ， 使 得 对 数据 的 描述 变 得 

棒 性 ,文献 [5] 提 出 奇异 值 分 解 与 小 波 变 换 相 结合 的 水 印 算法 。 方便 与 合理 ， 同 时 还 在 一 定 程度 上 抑制 外 界 变化 对 特征 提取 

先 将 图 像 进行 小 波 变换 , 对 低频 部 分 进行 SVD 后 将 水 印 内 容 。 造成 的 影响 ， 所 以 NMF 分 解 算法 相 较 于 传统 的 一 些 算法 而 

进行 嵌入 。 文 献 [6] 利 用 DCT 与 DWT 嵌入 鲁 棒 水 印 和 脆弱 水 言 ， 上 共有 实现 上 的 简便 性 、 分 解 形式 和 分 解 结果 上 的 可 解释 

印 实现 双 功 能 水 印 算法 ， 可 以 有 效 抵抗 压缩 、 噪 声 、 局 部 前 。 性， 以 及 占用 存储 空间 少 等 诸多 优点 。 

切 的 攻击 。 虽 然 文献 [5,6] 两 中 算法 均 满 足 水 印 算法 的 隐蔽 性 1.2 奇异 值 分 解 

和 重 棒 性 ， 但 算法 中 奇异 值 分 解 带 来 的 对 角 线 失真 以 及 虚 警 奇异 值 分 解 (SVD) 是 最 著名 ， 最 广泛 使 用 的 矩阵 分 解 方 

率 问 题 没有 得 到 良好 的 解决 。 文献 [7] 提 出 了 一 种 基于 奇异 值 。 法 。 经 过 奇异 值 分 解 后 ， 对 应 的 正 交 和 矩阵 表示 图 像 的 几何 结 

分 解 与 蜂 群 优化 的 鲁 棒 水 印 算法 ， 利 用 蜂 群 优 化 算法 来 选取 构 ， 奇 异 和 矩阵 表示 图 像 的 亮度 信息 。 因 该 分 解 算法 拥有 强 稳 

嵌入 强度 ， 自 适应 的 均衡 了 水 印 算法 的 鲁 棒 性 与 透明 性 , 但 定性， 被 广泛 地 应 用 于 图 像 处理 领 域 。 设 定数 字 图 像 矩 阵 

该 优化 算法 收敛 速度 慢 ， 寻 找 最 优 解 时 间 较 长 。 4sRxw ，R 为 实数 域 。 正 交 和 矩阵 seRww , Ven" SD RAE 
为 从 根本 上 解决 水 印 的 透明 性 与 鲁 棒 性 不 均衡 问题 , 温 E S eR 使 得 4 表示 为 

泉 等 人 四 提出 了 零 水 印 算法 ， 利 用 原始 载体 图 像 的 内 部 特征 A-UxSxVT Q) 

来 构造 零 水 印 。 在 保证 了 原始 载体 图 像 完整 性 同时 ， 又 有 效 其 中 UU 表示 左 奇异 值 正 交 矩阵 和 信 表示 右 奇异 值 正 交 和 矩阵 

地 解决 了 水 印 算法 鲁 棒 性 与 透明 性 的 矛盾 问题 。 曲 长 波 等 人 HAA E UU =1, V-V^. S 是 一 个 非 对 角 元 素 值 都 为 0 

中 提出 基于 Curvelet-DSVD 和 视觉 密码 的 强 鲁 棒 零 水 印 。 依 ”的 矩阵 满足 44=L12…7) 且 4>…>4>0 。 7 为 矩阵 4 的 秩 。 于 

据 视觉 密码 生成 的 共享 份 与 双 奇 异 值 分 解 得 到 的 特征 矩阵 异 EA 

或 生成 零 水 印 。 该 算法 在 抵抗 JPEG 压缩 、 噪 声 、 滤 波 、 剪 À 0 0 

切 方 面 都 表现 出 较 好 的 鲁 棒 性 ， 但 算法 复杂 且 在 抵抗 几何 攻 rio (4) 

击 方面 还 存在 很 大 的 提升 空间 。 在 虚 警 率 问 题 上 实验 结果 仍 I Bh 

不 理想 。 本 文 借鉴 文献 [10] 的 思路 提出 一 种 基于 分 块 的 NMF 4 是 由 分 解 而 唯一 确定 的 , 记 做 矩阵 4 的 奇异 值 ,分 解 式 
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区 域 不 重大 分 块 ， 
子 块 分 解 ， 得 到 


基数 入 


对 图 


像 进行 二 级 小 波 变 换 后 的 


时 与 经 过 Arnold 与 混沌 映 


E 负 特性 利用 NMF 对 不 
和 矩阵 


重 置 乱 加 密 水 印 图 


须 优化 算法 (BAS) 


算法 理论 基础 
非 负 和 矩阵 分 解 


放 比 例 的 参数 8。 该 算法 的 全 
印 算法 提取 时 出 现 的 高 虚 警 错误 和 稳健 性 不 强 问 题 


像 进行 异 或 运算 生成 零 水 印 。 利 ) 


UVxSxy7 记 做 4 的 


对 应 
像 结 


图 像 在 增强 奇异 


奇异 值 分 解 。 
值 分 解 后 ， 


的 关系 ， 不 同 的 图 


构 却 没有 任何 


IK 


自 适应 确定 增强 奇异 值 分 解 中 最 抗 攻 


搜索 算法 实现 自 适 


应 过 程 提升 鲁 棒 性 。 


i 势 在 于 解决 奇异 值 分 解 水 


E 
结合 


非 负 矩阵 分 解 (non-negative matrix factorization,NMF) 


[视觉 和 图 像 工 程 等 研究 领域 中 


行 水 
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字 图 
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欢迎 的 多 维 数据 处 理 , 


均 为 非 负数 约 束 条 件 下 的 矩阵 分 解 方法 。 算 法 描述 如 下 : 
找 非 负 基 和 矩阵 W e nn 和 
后 可 理解 为 : 原始 和 矩阵 


给 定 矩 阵 YeRe"，, 寻 
be nes, [fb vo wxn , SRN 


EHETE ERLA 


E 负 特征 矩 


向 量 是 对 左 矩 阵 亚 中 所 有 列 向 量 的 加 权 和 ， 而 权重 系数 


HIERE H 对 应 列 向 量 的 元 素 ， 
矩阵 或 特征 矩阵 。r 为 要 分 解 矩 
m 小 ， 即 满足 (m+n)r<mn , XXISEFH'EHGEXR REAR 
始 矩 阵 ， 就 可 以 实现 对 原始 和 矩阵 进行 降 维 
} NMF 矩阵 分 解 实 


要 比 n 和 


降 维 矩阵 。 


AE 


目标 
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减 。 
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函数 为 


贰 则 为 


NMF 与 PCA、 


Pe V IW HERE HX 
数 的 设 定 和 和 迭代 的 选择 。 定 义 如 下 : 


E 阵 的 乘积 


FA. ICA 等 常 | 


分 解 后 的 所 有 


故 称 W 为 基 和 矩阵 ， 
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分 量 均 为 非 负 值 并 且 实 现 非 线性 的 维 数 约 


这 种 非 负 性 的 限制 导致 了 相应 描述 在 一 定 程度 上 的 稀 下 C 


的 对 


H.rH. 


增强 奇异 值 分 解 
Wake FL 


ff A 大 小 为 MxM 


4=Ux(S) xVr,.0« «1 
: 矩阵 4 的 左 奇异 
增强 奇异 值 矩 分 解 是 将 数字 图 像 4 进行 
角 和 矩阵 Suo 做 8 次 肾 运 算 。 在 处 理 不 
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TOI 
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(BN-SVD) 03 
出 的 一 种 新 的 分 解 
。 其 作用 是 均衡 奇异 值 矩 阵 对 角 


Vus 的 正 交 和 矩阵 和 Suas 


图 像 与 奇异 值 向 量 构成 
有 相同 的 奇异 值 ， 但 两 幅 图 
因 这 一 缺点 导致 经 过 奇异 值 向 量 进 
的 图 像 中 提取 出 水 


像 可 以 


是 基于 奇异 值 分 解 的 理 
算法 ， 对 奇异 值 矩阵 进行 修改 引入 参 
线 方向 上 的 灰 度 。 设 数 
其 进行 增强 奇异 值 分 解 ， 则 存在 
的 对 角 和 矩阵 ,公式 成 立 如 下 : 


， 对 


(5) 


矩阵 是 UV， 右 奇 异 矩 阵 是 V7 。 
奇异 
司 的 数字 图 像 时 ， 
的 信息 对 参数 进行 合理 的 调整 ， 以 达 


佳 的 抗 攻击 效果 。 
值得 到 


AE 


存在 


放大 ， 降 低 
程度 上 提高 了 算法 
的 虚 警 错误 问题 上 


其 增强 奇异 值 分 解 的 优势 在 于 : 图 像 
图 像 矩 阵 在 受到 攻击 时 的 敏感 性 ， 

的 鲁 棒 性 。 而 且 在 解决 奇异 值 分 解 中 

， 增 强 奇 异 值 分 解 对 奇异 值 向 量 做 特 
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像 的 特征 ,这 样 的 处 理 方式 使 虚 警 问题 得 到 更 好 的 解决 。 
自 适 应 天 牛 须 优化 算法 


天 和 牛 须 (beetle antennae search, BAS) 算 法 
虫 搜索 行为 启发 , Æ 2017 年 提出 了 一 种 基于 甲虫 砚 食 原 
优化 算法 。 与 粒子 群 ， 遗 传 算法 相似 ， 无 须知 道 
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这 一 简单 原理 天 


牛 就 可 以 有 效 找到 食物 。 


算法 原理 如 下 : 
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录用 定稿 


a) 初 始 化 随机 天 牛 位 置 气味 强度 x 和 天 牛 须 朝 向 的 随机 
向 量 b 并 对 此 归 一 化 ， 设 置 最 大 迭代 次 数 MAX o 
> rands(k,l 
Ti m 
其 中 : rands0 为 随机 函数 ; 上 表示 空间 维度 。 

b) 天 牛 个 体 两 须 采 集 自身 附近 两 点 的 气味 强度 Xw X, 
并 朝向 食物 的 坐标 。 


X,-X'«d,xb]h 
人 (7) 
X,- X'-d,xBfo 
其 中 :表示 在 第 1 次 迭代 时 ， 天 和牛 左右 须 位 置 坐标 为 Xo X, ; 
天 牛 质 心 坐标 X'; 两 须 之 间 的 距离 do o 
c) 和 迭代 更 新 天 牛 位 置 的 气味 强度 X. 


X 2 X! 2à xbxsign(f(X,, )- fx, ) (8) 


EP: o 表示 在 1 次 迭代 时 的 步 长 因子 ;sign0 为 符号 函数 


d) 由 于 参数 Ap ， 取 值 范围 属于 [0,1]， 即 在 本 实验 中 气味 
强度 X 需要 进行 分 段 约束 。 
9 XU T 
X -40 X" «0 (9) 
X=X" else 


e) 天 牛 个 体 两 须 附近 两 点 的 气味 强度 坐标 ( Xs , X, ) 利用 
式 (9) 确 定 当前 天 牛 空 间 位 置 气味 强度 和 气味 强度 测定 函数 
f 〔 即 适应 度 函数 ) 来 求 出 天 牛 个 体位 置 SX): 


f(X)=f (x) (10) 
有 通过 测定 函数 f 寻 到 天 牛 个 体位 置 两 须 中 最 优 的 气味 
强度 X o 
[AX] 9 max(X) (11) 
外 测定 天 牛 位 置 最 强 气 味 强 度 Xbest 并 探寻 到 最 优 测定 
函数 jbest 奔 向 该 位 置 。 
f = max(fbest) 
s = max(Xbest) (12) 
hb) 重 复 执行 步 又 b)-f£), 3SRUUIEOIE SEHR IS DUAE, ERR 


和 迭代 结束 时 ， 判 断气 味 ， 如 果 当 前 气味 强度 优 于 前 一 次 迭代 
的 好 则 执行 g)， 当 执行 最 大 迭代 次 数 MAX 则 结束 整个 过 程 ， 
得 到 最 优 结果 xX BEA A 。 
2 ”数字 水 印 算法 实现 
2.1 水 印 图 像 的 预 处 理 

Arnold 变换 5 俗称 猫 脸 映射 是 一 种 从 规则 位 置 到 随即 
立 置 的 映射 ,是 一 种 传统 的 混沌 系统 ,对 一 副 大 小 为 MxM 的 
图 像 ，Arnold 置 乱 定 义 如 下 : 


B 7 l i 四 (RED (3) 


x, y e(0,L--N-1) 


其 中 前 后 像素 点 的 位 置 分 别 由 (xsy) 和 (x,y) 表示 。w%yN 为 正 
整数 ，WN 代表 图 像 矩 阵 的 阶 数 ， mod 为 取 余 运算 。 其 变换 的 
实质 是 改变 像素 位 置 ， 来 减弱 相 邻 像 素 之 间 的 关联 性 ， 提 高 
了 图 像 安 全 性 。 拥 有 周期 性 ， 经 过 TIR e Fe T M Hi 
原始 图 像 。 

Logistic 映射 49 具有 遍历 性 高 、 
感 型 等 优点 。 其 定义 如 下 : 
x 7 tx(17x ) x; € (0,1) u e[0,4] (14) 


伪 随 机 性 强 对 初 值 高 敏 


xs(0D IM, Logistic 映射 工作 处 于 混沌 状态 ， 也 就 是 说 ， 
初始 条 件 % 在 Logistic 映射 作用 下 产生 的 序列 是 非 周期 性 


ay LK m 


肖 振 久 ， 等 : NMF 和 增强 奇异 值 分 解 的 自 适 应 零 水 印 算法 


中 :4 为 Logistic 混沌 系统 的 控制 参数 , * 为 映射 的 混沌 序列 。 
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的 ， 不 收敛 的 ， 而 在 此 范围 之 外 ， 生 成 的 序列 必 将 收敛 于 某 
一 个 特定 的 值 ， 当 3.569456<u<4 时 ， 特 别 是 比较 靠近 4 时 ， 
迁 代 生成 的 值 是 处 于 一 种 伪 随 机 分 布 状态 即 混 沌 状态 ， 而 在 
取 其 他 值 时 ， 在 经 过 一 定 次 数 的 迭代 之 后 ， 生 成 的 值 将 收敛 
到 一 特定 的 数值 。 
2.2 零 水 印 构 造 
本 文采 用 离散 小 波 变 换算 法 作 基 底 ， 对 原 图 像 进行 二 级 
离散 小 波 变 化 。 因 原始 图 像 能 量 高 度 集中 于 图 中 位 于 低频 频 
段 的 LL2 区 域 ， 而 细节 信息 则 主要 分 布 于 高 频频 段 。 由 于 
低频 系数 稳定 性 高 ， 所 以 选择 在 其 中 供 入 水 印信 息 ， 保 证 水 
印 的 强 鲁 棒 性 且 抗 攻击 能 力 得 到 大 大 提升 。 


原始 灰 度 图 A Eos r 
= J 4 Hj 
像 I 近 子 图 LL2 diuinum 水 印 图 像 W 

| 
JE ELES) EN K—— K 
兴隆 分 一 一 | za ! 
AEREE BN-SVD 计 算 最 大 奇 奇 
特征 矩阵 Hi 一 一 一 一 一 3 异 值 均值 e k 
Logisitic 映 射 


加 密 后 水 印 W 


— 特征 向 量 F 上 一 


零 水 印 图 像 M 


图 1 零 水印 构 建 图 
Fig. 1 Zero watermark construction diagram 

a) 大 小 MXM 原始 灰 度 图 像 上 /进行 2 级 DWT; 对 2 级 

DWT 变换 后 的 得 到 的 低频 逼近 子 带 LL2 HIT nxn 分 块 ， 生 

BUR SGRTHR A [i-is DL], EXE A HEITA r 的 NMF 


Axnxn 
分 解 ， 


Hrzn , 即 


A, «W,H, (15) 
其 中 : W, 和 1; 分 别 是 mxr 0 rxn BE fA yE RE 
b) 对 特征 矩阵 H: 进行 增强 奇异 值 分 解 ， 即 


H, -U,x (S) V7 (0 f «1) (16) 


其 中 : Ui,V; 为 正 交 矩阵 5; 为 对 角 矩 阵 。 其 参数 5 用 上 述 1.4 
天 牛 须 搜 索 找到 最 优 。 

c) BOSE fü B Ee S; 中 第 一 个 奇异 值 共有 (M/2n)x(M/2n) 个 ， 
记 作 m L2, (MJ2n)x(MJ2n)) 


由) 计算 (Mj2n)x (M[2n) 个 最 大 奇异 值 的 平均 值 meant, ) 。 
依据 六 和 meann) 的 大 小 关系 构成 一 个 特征 向 量 F 。 


1 1;)mean(n,;) 
H- P e (17) 
e) 将 原始 二 值 水 印 图 像 W 作 Ki 次 Arnold 置 乱 ， 得 到 置 


乱 后 水 印 图 像 W,。 然 后 对 水 印 图 像 W, 进行 Logistic 映射 ， 经 
过 二 次 加 密 后 水 印 图像 W, 和 密 钥 K 。 


f) 将 特征 向 量 F 与 置 乱 后 的 水 印 图 像 WW 进行 异 或 运算 ， 
生成 零 水 印 图 像 4M : 
M - XOR(F.W,) (18) 


2.8 零 水 印 检测 
水 印 的 检测 即 为 水 印 构造 的 逆 过 程 , 如 图 2 所 示 。 
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0.5070，0.1382， 达 到 强 鲁 棒 性 效果 。 分 别 对 三 幅 灰 色 图 像 
进行 了 几何 攻击 和 非 几何 攻击 ,图 4 为 部 分 表 1 中 Lean 图 像 
受到 不 同类 型 攻击 后 的 实验 效果 图 ， 表 1 列 出 了 待 测 图 像 在 
受到 各 种 不 同 攻击 后 提取 出 32x32 大 小 水 印 图 像 的 NC 值 
从 表 中 数据 可 以 得 出 ， 本 算法 可 以 有 效 抵抗 各 种 攻击 ，NC 
值 的 平均 值 高 于 0.9984， 随 着 攻击 强度 的 增加 和 混合 攻击 ， 
水 印 的 NC 值 受到 一 定 的 影响 。 其 中 ， 当 旋转 攻击 强度 为 10 
度 时 ， 提 取出 的 水 印 的 NC 值 最 低 为 0.9988。 由 此 可 见 本 文 


录用 定稿 
原始 灰 度 图 2 级 低频 通 
[TUN 近 子 图 LL2 | 让 分 块 处 理 - 
— I 水 印 图 像 W 
非 负 短 阵 分 eg 
AR (NMF) FIRA; 
| 逆 Anor1ld 变 换 
BN-SVD [agya se 
EE 计算 最 大 奇异 
TEAERBIEH mias de, 
i nn 道 Logisitic 映 射 
Fi=1l Fi=0 "n 
IRE Jr a JAK ERW 
一 一 特征 向 量 F " We 
XOR 
零 水 印 图 像 M 


图 2 零 水 印 检测 图 

Fig.2 Zero watermark detection diagram 

选取 原始 灰 度 图 像 T' ， 大 小 为 MxM 。 零 水 印 检 测 步 又 
如 下 : 

a) 对 原始 灰 度 图 像 I' 进 行 2 R DWT, 对 得 到 的 低频 逼近 


子 带 LL2 进行 axn 分 块 ,生成 不 重 肥 子 块 , (i=12.… 


并 对 A; HEITER r 的 NMF AE, Horn , 即 
A xW; H; (19) 
b) 对 特征 矩阵 H; 进行 增强 奇异 值 分 解 (BN-SVD)， 得 到 
对 角 和 矩阵 S; 。 
cX AERE S; 中 第 一 个 奇异 值 共 (M/2n)x (MJ2n) 个 记 
作 : m(G-L2-.(M/2n)x(MJ2n)) 。 
d) 求 取 平 均值 meann; ) 。 根 据 n 和 meanfyw ) 的 大 小 关系 得 
到 特征 向 量 F'。 
6) 将 零 水 印 图 像 M 与 人 中 得 到 的 特征 向 量 ' 进行 异 或 
运算 ， 得 到 加 密 之 乱 后 的 水 印 图 像 W o 
W, 2XOR(M,F") (20) 
名 最 后 对 水 印 图 像 更 进行 道 Logistic BU] AU K, 次 Arnold 
逆 置 乱 得 到 原始 水 印 图 像 灰 。 


3 ”仿真 实验 结果 与 分 析 


本 实验 主要 在 matlab2016a 实验 平台 下 仿真 ， 实 验 中 选 
取 图 3(a)~(c)Lena,Baboon,Plane 大 小 均 为 512x512 像素 的 标 
准 灰 度 图 像 作 为 测试 图 像 ， 水 印 图 像 选取 图 (d) 大 小 为 32x32 
像素 的 二 值 图 像 “ 辽 宁 工 大 ”。 设 置 天 牛 须 优化 迭代 次 数 为 
20。 通 常 水 印 算法 的 评价 指标 是 不 可 见 行 和 和 鲁 棒 性 ， 但 在 本 
文中 ， 对 于 零 水 印 算法 ， 水 印信 息 没 鹏 入 到 载体 图 像 中 ， 
即 需要 和 鲁 棒 性 作为 评价 指标 ,仿真 实验 采用 归 一 化 相关 (NC) 
函数 来 评价 提取 的 水 印 与 原水 印 的 相似 度 , 即 

YY hol) 


NC(x1,x2)= d Q1) 


M N MN 
Y uto jo I x jp 
i=] j=l i=] j=l 


其 中 : 如 表示 为 初始 水 印 ， 志 表示 为 提取 的 水 印 
3.1 和 鲁 棒 性 实验 

为 验证 本 文 水 印 算法 具有 较 好 的 鲁 棒 性 ， 本 文采 用 NC 
作为 算法 鲁 棒 性 的 评价 标准 。 针 对 32X32 大 小 的 水 印 图 像 ， 
本 实验 利用 天 牛 须 优化 算法 在 第 16 次 实验 中 得 到 Lena, 
Baboon 和 Plane 三 幅 灰 色 图 像 的 最 优 参 数 p. 值 为 0.4289， 


， 
Axnxn 


E 
, 


算法 可 以 有 效 抵抗 旋转 攻击 。 对 于 前 切 攻击 ， 前 切面 积 增 大 
提取 出 的 水 印 NC 值 在 0.9984 以 上 ; 噪声、 滤波、 缩放 、 压 
缩 和 混合 攻击 在 强度 增加 时 ， 提 取出 的 水 印 NC 值 都 接近 于 1。 


(b)Baboon 


[ 


辽宁 
IX 


(c)Plane 


(d)32x32 水 印 图 像 
图 3 原始 灰 度 图 像 与 水 印 图 像 
Fig.3 Original grayscale image and watermark image 
分 析 以 上 的 结果 可 得 ， 通 过 NMF 和 BN-SVD 结合 的 算 
法 ， 并 利用 天 牛 须 优化 算法 智能 搜索 最 佳 的 参数 8， 增强 了 
该 水 印 算法 的 鲁 棒 性 。 
表 1 待 测 图 像 攻击 后 提取 出 32X32 水 印 的 NC 值 


Table 1 Extracting the NC value of the 32x32 watermark after 
attacking the image to be detected 
攻击 强度 Lena Baboon Plane 
言 斯 噪声 0.01 0.9999 0.9998 1.0000 
abd 0.02 0.9999 0.9999 1.0000 
3x3 1.0000 0.9999 1.0000 
中 值 滤波 
5x5 0.9999 0.9999 1.0000 
高 斯 滤波 0.5 1.0000 0.9999 1.0000 
20 1.0000 1.0000 1.0000 
压缩 攻击 
10 1.0000 0.9999 1.0000 
缩放 0.5 1.0000 0.9999 1.0000 
高 斯 滤波 加 中 
混合 0. 0.9988 0.9990 0.9990 
q 899] 1/4 
左上 角 1/16 1.0000 0.9999 1.0000 
剪 切 攻击 " 
左上 角 v4 0.9990 0.9984 0.9992 
椒盐 噪声 0.01 1.0000 0.9999 1.0000 
0.02 1.0000 0.9999 1.0000 
向 左 10° 0.9988 0.9990 0.9990 
旋转 攻击 
向 右 10° 0.9987 0.9990 0.9990 
3.2 ”安全 性 分 析 
在 整个 算法 中 ， 安 全 性 是 不 容 忽视 的 。 本 文 算法 对 水 印 


图 像 进行 Arnold 变换 和 Logistic 映射 ,实现 水 印 图 像 特征 的 二 
次 加 密 , 加 密 密 匙 分 别 为 K，K, ,这 样 即 使 零 水 印信 息 泄露 ， 
在 不 确定 密 匙 的 情况 下 ， 是 无 法 恢复 出 水 印信 息 的 ， 说 明 本 
算法 的 安全 性 较 高 。 
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(a) 椒 盐 噪 声 (0.02) 
(a) Salt and pepper noise (0.02) 
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NT 


LA 


(b) 高 斯 噪声 (0.02) 
(b)Gaussian noise(0.02) 


ALT 辽宁 
工大 工大 
(Ac E fT 1/4 (d) 高 斯 滤波 加 中 心 剪 切 4 


(c)Shear upper left corner(1/4) 


I MEE aa tm 


(a) Ug EK; 错误 提取 的 水 印 


(a) Extracted watermark of 


辽宁 
工大 
(e)JEPG 压缩 (20) 
(e)JEPG compression(20) 
图 4 受 攻击 后 的 实验 图 像 


(d)Gaussian noise and Central 


shear(1/4) 


辽宁 
工大 


(Oh (EISE GS) 


(f)Median filtering(5x5) 


Fig.4 Attacked experimental images 


only the error key K2 


ione 


A3 nc 


ORL K, 和 K 均 错 误 所 提取 的 水 印 


muu 5 
L] 
"E 


FEL 


(b) 仅 密 匙 


(b)Extracted watermark of 


由 3 


i MEE 
uriy: E 
"IT Ee 


K, 错误 提取 的 水 印 


only the error key K1 


(c)Extracted watermark of only the error key K1 and error key K2 


图 5 不 同 密 钥 提取 的 错误 


吴 水 印 


ESI 


[spun 


Fig.5 False watermark information extracted from different keys 


图 


5 可 知 在 密 匙 错误 的 情况 下 提取 出 的 水 印 情 况 ， 图 


5 中 , (a) 为 在 密 匙 Ki 正确 Ko 错误 的 情况 下 提取 出 的 水 印信 息 ; 


(b) 为 在 密 匙 K, 正确 Ki, 错误 的 情况 下 提取 出 的 水 印信 息 ; (c) 


HEK RU 


EX] ut n] 


b K, K, 都 错误 的 情况 下 提取 出 的 水 印信 息 。 从 图 5 


中 可 以 得 到 ， 在 密 古 错误 的 情况 下 是 无 法 提取 出 水 印信 息 ， 


以 表明 本 文 算法 的 安全 性 较 高 。 


3.3 虚 警 率 问题 实验 测试 


为 


误 问 题 


了 验证 本 文 是 否 能 够 真实 有 效 的 解决 SVD 的 虚 警 错 
， 将 本 文 算法 与 文献 [17] 和 传统 SVD 算法 进行 对 比 实 


验 , 图 6 中 , (a) 为 原始 灰 度 图 像 ，(b) 为 伪造 的 水 印 图 像 ，(c) 


为 本 文 


算法 提取 出 的 水 印 。 实 验 表明 , 本 算法 引入 参数 8 后 ， 


便 存 在 


图 像 与 奇异 值 向 量 一 一 对 应 的 关系 ， 在 不 确定 参数 8 


值 的 情况 下 不 能 正确 的 提取 出 水 印信 息 。 


RRE 
beth 


RRE 
beth 


RRE 
beth 


(a) 原 始 图 像 (b) 伪 造 水 印 图 像 (0) 本 文 算法 
(a) original image (b)forgery watermark 提取 出 的 水 印 
image (c)the algorithm 


in this paper 


图 6 虚 警 错误 仿真 实验 结果 


Fig.6 Simulation results of false positive errors 


表 2 中 三 种 算法 都 采用 离散 小 波 低频 分 块 ， 区 别 在 于 对 
和 矩阵 使 用 了 不 同 的 分 解 方法 。 其 中 文献 [17] 提 出 了 Schur 算 
法 。 由 表 2 可 知 实验 后 Schur 算法 与 本 文 算 法 提取 的 水 印 图 
像 与 原始 水 印 图 像 的 NC 值 都 在 0.5 以 下 ，SVD 算法 虚 警 问 
题 明 显 居 高 ， 但 本 文 算法 NC 值 在 0.4 以 下 ， 说 明 本 文 算法 
处 理 虚 警 问题 能 力 强 。 
表 2 ”本文 算法 与 文献 [17]NC 值 对 比 
Table 2 The NC value of the algorithm is compared with the literature 
[17] 
图 像 SVD 算法 文献 [17] 本 文 算法 
Lena 0.9963 0.4999 0.3558 
Bridge 0.9896 0.4875 0.3146 
Plane 0.9942 0.4968 0.3348 


3.4 本 文 算法 与 文献 [10,18,19] 的 对 比 实验 


为 
级 低频 
对 特征 


了 进一步 验证 本 算法 的 优势 ， 本 算法 对 DWT 后 的 二 
区 域 分 块 利用 NMF 算法 分 解 后 得 到 的 特征 矩阵 ， 并 
矩阵 进行 增强 奇异 值 分 解 ， 通 过 分 解 后 均衡 化 对 角 线 


方向 上 
自 适 应 


的 灰 度 ， 使 水 印 拥有 强 鲁 棒 性 ， 利 用 天 牛 须 优化 算法 
确定 参数 8 ， 使 水 印 的 鲁 棒 性 达到 最 佳 。 本 文选 取 


Lena 图 像 作为 原始 灰 度 图 像 , 与 文献 [10,18,19] 算 法 得 到 的 实 


验 结果 
EAR 
过 细胞 


进行 了 对 比 ; 文献 [10] 采 用 的 是 离散 小 波 低频 分 块 后 
余弦 变换 再 通过 增强 奇异 值 分 解构 成 特征 向 量 ， 并 通 
神经 网 络 确定 可 逆 逻 辑 运 算 构成 零 水 印 ， 文献 [18] 是 
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在 时 域 上 使 用 非 均匀 NURP 构造 特征 矩阵 ， 并 通过 SURF 来 的 问题 ， 此 外 在 原始 灰 度 图 像 受到 几何 攻击 与 非 几何 攻击 下 
对 图 像 进行 几何 校正 ， 文 献 [19] 对 像素 近似 无 损 的 安全 区 域 。 所 提取 出 的 水 印 NC 值 均 可 达到 98% 以 上 ， 增 强 了 算法 的 鲁 
进行 二 级 元 余 离 散 小 波 变 换 的 低频 区 域 分 块 ,再 SVD 分 解构 — 棒 性 。 然 而 本 文 算法 在 压缩 攻击 和 缩放 攻击 方面 仍 需 进 一 步 
造 特征 和 矩阵， 而 本 文 算法 采用 的 是 二 级 离散 小 波 低频 分 块 后 。 探索 。 
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